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ABSTRACT 


The  short  time  response  of  a  cylindrical,  elastic  segment 
to  shock-wave  engulfment  is  treated.  The  exact  equations  of  the 
linear  theory  of  elasticity  are  used  instead  of  the  usual  integrated 
equations  of  ring  theory.  Using  a  relatively  new  mathematical 
technique,  described  in  the  report,  we  easily  identify  the  focusing 
effects  due  to  the  curved  geometry.  This  work  explains  some 
anomalous  results  obtained  while  subjecting  curved  aluminum 
specimens  to  loads  produced  by  shock  waves. 
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I.  INTRODUCTION 

Recently,  a  facility  was  developed  to  measure  the  short  time  effects 
of  the  response  of  materials  to  a  step-function  pressure  pulse  [  l] .  This 
has  proven  efficacious  in  the  measurement  and  prediction  of  the  dispersive 
properties  of  materials  [lO].  Loading  is  applied  to  the  specimen*  by  means 
of  a  shock  tube  and  provides  a  clean  step  function  of  up  to  70  psi  with  a  rise 
time  of  less  than  15  nsec.  The  back-face  velocity  of  the  specimen  is  measured 
by  a  capacitance  gauge,  and,  in  the  usual  case,  the  one-dimensional  theories 
used  are  good  for  a  few  reverberations.  A  flat  aluminum  specimen  is  used 
for  calibration  since  it  is  quite  dispersion  free.  A  typical  oscilloscope  trace 
of  the  back-face  velocity  is  shown  in  Fig.  1. 

In  order  to  ov'-mu  this  technique  to  geometries  other  than  flat  specimens, 
recently,  some  cylindrical  segments  ware  tested.  The  results  for  an  aluminum 
specimen  {Fig.  Z)  show,  instead  of  the  usual  "staircase"  response  (Fig.  1), 
a  jump  followed  by  a  rise  and  then  a  downward  slope.  It  was  easily  ascer¬ 
tained  that  this  phenomenon  was  not  due  to  malfunctions  of  the  test  apparatus. 
Thus,  this  effect  must  come  from  the  curved  geometry. 

It  is  the  purpose  of  this  woi'k  to  explain  these  phenomena.  Two  major 
assumptions  were  used  in  order  to  simplify  the  problem.  The  first  was  to 
idealize  the  ring  segment  to  that  of  a  full  ring.  This  was  a  valid  approxima¬ 
tion  since  only  short  times  were  of  interest,  and  the  experiment  was  normally 

Although  the  specimen  can  vary  in  size,  the  usual  sample  is  a  flat  specimen 
1 . 25  in.  X  1 . 25  in.  X  0.  25  in. 
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Fig.  2.  Curved  Specimen  Trace 
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considered  finished  when  effects  from  the  boundaries  of  the  specimen  arrived 
at  the  point  where  the  back-face  velocity  was  measured,  in  this  case,  at 
0^0  (Fig.  3).  The  second  assumption  was  that  the  effects  of  shock  diffrac¬ 
tion  were  negligible,  and  that  the  shock  wave  acted  as  a  pressure  pulse, 
unchanged  in  pressure  as  it  engulfed  th'  specimen.  This  assumption  was 
made  on  the  basis  of  preliminary  calculations  and  the  fact  that  the  theory 
compared  well  with  experiments. 

The  results  of  the  present  analysis  show  that  the  first  ramp  is  due  to 
successive  P  wave  arrivals  as  the  engulfment  ensued  and  to  the  focusing 
effects  of  curved  media.  The  downward  ramp  is  shown  to  be  due  to  S  wave 
arrivals.  Rather  than  the  usual  modal  approach,  we  proceed  with  the  tech¬ 
nique  used  by  Friedlandcr  [4],  in  which  the  solution  is  taken  as  a  wave  sum, 
each  term  of  which  is  identified  with  one  sheet  of  an  infinite-sheeted  Riemann 
surface  over  the  propagation  coordinate  (circumferential  0  here). 

The  problem  of  the  engulfment  of  a  ring  has  been  treated  by  several 
authors  of  which  Bleich  and  Mindlin  [2]  and  Payton  [3]  are  the  most  germane 
to  the  present  work.  These  authors  are  primarily  interested  in  the  structuial 
response  of  a  ring  rather  than  the  wave  propagation  through  the  thickness  as 
in  the  present  work.  The  equations  of  the  linear  theory  of  elasticity  are  the 
equations  pertinent  to  the  present  problem  and  not  the  integrated  ring 
equations  used  in  references  [2]  and  [3j. 
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II .  ANALYSIS 


As  was  mentioned  in  the  Introduction,  a  cylindrical  specimen  is  hit  by 
a  shock  wave  that  is,  to  a  close  approximation,  a  step-function  pressure 
input.  The  problem  considered  in  the  present  work  is  the  calculation  of  the 
velocity  of  the  back  surface  directly  behind  the  specimen  at  0  =  0.  The  times 
of  interest  are  those  necessary  for  the  first  few  reverberations.  Accordingly, 
a  solution  valid  for  small  times  is  sought. 

Consider  an  elastic  ring  of  inner  radius  r^  and  outer  radius  r^.  It 
has  a  flat  pressure  pulse  of  strength  pQ  impinging  on  it  with  velocity  V 
(Fig.  3).  The  elastic  parameters  are  X.,  p,  and  p,  and  the  problem  is  treated 
as  that  of  plane  strain.  Time  is  denoted  by  t,  and  polar  coordinates  by  r 
and  0.  Setting 


r 


(1) 


where  a  and  (3  are  the  dilatational  and  equivoluminal  speeds  given  by 
2  2 

a  -  (X.  +  2p)/p  and  (3  '  =  p/p,  respectively,  one  gets  the  equations  of  elasto- 
dynamics 


vWtt  vW*,tt 


rjU  =  <S>,r  +  r'1^,  Q  r<v  =  r"V 


1 


e  r’  r 


(2) 

(3) 
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rl^r  =  (t,2  -  2)V2£  x  2(0,  rr  -f  r'^,xQ  -  r‘2^0) 


V  “  2(r  r0  "  0)  ~  ^  rr  *  r  r  +  r 


4-  r'^tb 


-2 


96 


rl(<rr+,r0)  2  2 

-J— 7; - “=  2{ti^  -  1)V> 


(4) 


where  V2  =  92/3r2  +  r~*3/3r  +  r_23/302,  u  and  v  are  the  radial  and  circum¬ 
ferential  displacements,  respectively,  and  c r  ,  <Tg,  and  (r^g  are  the  usual  stress 
components.  The  initial  and  boundary  conditions  are 

0=^.t  =  ‘j>  =  4J,t=O  y<r<l  t  =  0  YsJl  (5) 

cr  —  -p  H(vt  -  r  cos  G)  0-0  =  0  r=i  t>0 

r  1  o  r0 

o’r  =  o’j.q  =  °  r  =  y  t  >  0  (6) 

t 

where  v  =  V la. 

In  theory,  one  can  solve  this  problem  by  expanding  all  quantities  in  a 
Fourier  series  in  0,  taking  a  Laplace  transform  in  time,  solving  the  resulting 
ordinary  differential  equation,  and  inverting  the  transform  by  residue  theory. 
Such  a  process,  however,  is  not  well  suited  to  short  times  or  to  discontinuous 
-olutions  [2] .  A  method  that  appears  to  be  particularly  well  suited  to  pro¬ 
blems  of  this  type  is  that  based  on  the  work  of  Friedlander  [4].  In  this  method, 


solutions  are  sought  of  the  form 


0,6(r,  0,  t)  =  ^  r,  0  +  2n~,  t)  (7) 

n=-w 

Instead  of  summing  the  "modes"  of  the  solution,  the  form  of  equation  (?) 
replaces  the  ring  by  a  fictitious  series  of  Riemann  sheets  spiralling  over  each 
other  to  provide,  an  analytical  continuation  of  0  to  ±®.  Longer  times  mean 
that  more  Riemann  sheets  must  be  looked  at.  However,  until  the  disturbance 
reaches  9  =  ±ir,  only  the  first  sheet,  n  =  0,  needs  to  be  considered.  In  the 
present  analysis,  this  is  ideal,  since  only  small  times  are  of  interest.  Either 
of  the  preceding  techniques,  when  carried  to  completion,  yields  the  exact  solu¬ 
tion.  The  method  of  Friedlander  is  used  here  to  provide  the  information 
desired  in  the  easiest  manner. 

Supplanting  the  Fourier  series  by  the  successive  infinite  sets  of  Riemann 
sheets,  we  extend  consideration  of  0  from  [-ir,  tt]  to  (-co,  ®).  Accordingly,  an 
exponential  1*  ourier  transform  is  indicated.  If  F(0)  is  a  function  and  F(v)  its 
transform,  then  the  pair  can  be  defined  as 

F(v)  =  /*'°F(0)  e1V0  d9 

*'-00 


F(0)  = 


2tt 


-iv0 

e 


dv 


(8) 


For  a  more  complete  explanation  of  this  method,  see  references  [  3 J  -  [6). 


Performing  this  operation,  and  at  the  same  time  taking  a  Laplace  transform 
in  time  on  equations  (2)- (4),  one  has 

Arrf'1i,I-cV2tszi?=o 


7  .  -It  ,2-2^22  =  A 

4»,  rr  T  r  r  '  ^  r  +  n  s  )«{»  =  0 


=  =  .  -1= 
r ,u  =6,  -  ivr  ib 

1  r  r 


r.v  =  -  ivr  *6  -  tb, 

1  •  -  r 


2  Jr  /  2  2  v2\  =  ^’r  .  / 

■i  s  +27j,-2  — ,21V(- 


ij/_  r 

2  r  i 


2  ^rO 
1  P 


=  2i"0>  -  41)  -  f 


22,  2v2\=  ^'r 

s  +  — =■•  +  2 - 

r  Z! 


2CTr+<r9)  ,,  2  2= 

- - - =2(1  -  1)»  0 


In  equations  (9)-(l  1 )-  s  denotes  the  Laplace  transform  variable  and  a  bar 


denotes  a  Laplace  transform. 

The  solution  to  equation  (9)  is 


<p  =  AIv(sr)  +  BKv(sr) 


ib  =  CIv(qsr)  +  DKv(qsr) 
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where  and  are  modified  functions  of  the  first  and  second  kinds,  and  A, 
B,  C,  and  D  are  constants  of  integration.  Equations  (10)  and  (11)  yield 


=  r  /  ,1  r  Iv(Tlsr)  Kvfrsr>] 

rju  =  s  [AIv(sr)  +  BK,(sr)J  -  ivqs  [C  +  D  -  — -j 

=  ivs  |A  —s—~  +  B  — ^ — J  -  qs  |ci'  fasr)  +  DK' (tisr)J 


r2  —  =  s2[AF1(sr)  +  BK2(sr)]  +  2ivr12s2[  CG1  (qsr)  +  DG2(qsr)] 


r2^  =  2ivs^[  AG^(sr)  -5-  BG2(sr)]  -  ti2s2[  CHj  (qsr)  +  DH2(qsr)]  (13) 


where 


Fj(x) 


2  2v^\  *v(x) 

q  +  — 9“jlv(x)  "  2 


x 


X 


F2(x) 


K'  (x) 
v '  ' 


q2  +  2^-)Kv(x)-  2  x 


M*)  KM 


Gi(x)=  —  - 

X 


X 


K  (x)  K'(x) 

G2<*» - 2 - — 


X 
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‘l<*>  =  ( 


2v2\  Mx> 

1  +  —5-)  I  (x)  -  2  — - 

2  /  v'  x 


H2(x) 


=  (**^w 


K„!x) 


The  boundary  conditions,  in  general  terms,  are 


=  -f(s,v) 

ffr0  =0  (at  r  =  1) 

=  ?r0  =  0 

(15) 

Thus,  the  boundary  conditions  yield 


f(v,  s)r' 


FjfsJA  +  F2(s)B  +  21vt12[G1(tis)C  +  G2(tis)D]  = - ^ 


ZivCG^sJA  +  G2(s)B]  -  112[H1(ns)C  +  H2(r,s)D]  =  0 
FjtysJA  +  F2(ys)B  +  2iwr|2[  G1  (t1Ys)C  +  G2(r)ys)D]=  0 

ZivfG^ysJA  +  G2(ys)B]  -  r)2[ H1  (riysJC  +  H2(r|Ys)D]  =  0  (16) 

Ideally,  equation  (13)  can  be  solved  for  A,  B,  C,  and  D,  and  then  explicit 
expressions  can  be  obtained  for  the  double  transform  of  all  pertinent  quantities. 
These  transforms  then  lead  to  the  formal  solution  (inversion  integrals)  in  the 
usual  way.  Since  the  present  work  deals  only  with  short  times  and  a  different 
method  is  used,  the  resulting  expressions  are  not  given. 
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In  order  to  apply  this  present  method,  the  input  or  pressure  function 
must  be  expanded  in  a  series  of  the  form  given  by  equation  (7).  This  can  be 
done  directly  [6]  by  applying  the  Poisson  summation  formula 


(17) 


The  applied  stress  at  r  =  1  is  given  by  the  first  equation  of  equation  (6).  Thus, 


cr 

T* 


-{r/v)  cos  0 
e 


(18) 


Following  reference  [  6]  and  using  equations  (17)  and  (16),  one  gets 


f(v,  S) 


(19) 


It  will  be  seen  later  that  the  absolute  value  sign  on  I  j  can  be  dropped. 

If  one  applies  the  Laplace  transform  theory,  the  behavior  of  a  function 
for  small  values  of  t  is  associated  with  large  values  of  s  [7].  However,  for 
the  present  problem,  in  order  to  invert  the  Fourier  transform,  an  integration 
with  respect  to  v  must  be  made  over  the  entire  range  v«(-co,  co).  One  is  thus 
led  to  consider  asymptotic  expansions  of  I(/  and  K  uniformly  valid  in  v 
foi  large  values  of  s.  Such  asympotic  expansions  have  been  investigated  by 
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Olver  [8,9]  and  are  [JO]: 


-^rh(s)(Uj  cosh  £  t  Vj  sinh  e)  +  (U^  cosh  6  +  sinh&)  =  0 


hfe)  vi  -  -TTT  h‘w)v2  =  0 

T1  V  S 


~^r  h(vs)U  .  +  rr11  .  u‘  =  0 
_y2g4  'Y  i  h(riYs)  2 


where,  for  convenience,  one  has  set 


,  A  cg(YS)  .  B  fir  e-sM  u  =  08(nvs)  +  D  fir  e-gMs) 

1  /i;  ',2  2  ys  V2 


v  =^_eS(Vs)  +  Bx/|<,-g(vs) 

1  ^ 


v  =_C_eg(^s)  .  D  fl  e-gtnvs) 

2  v/2tt  V2 


and  where 


£  =  g(s )  -  g(ys) 

6  =  g(ns)  -  gOvYs) 


In  order  to  assess  magnitudes  of  terms,  it  is  necessary  for  one  to  examine 
the  properties  of  £  and  6  before  solving  equation  (22).  It  suffices  for  one  to 
consider  e  alone. 


,  \  ,2  2,1/2  ,  2  ,  v2  2,1/2  ,  „  vr(v2  +  s2)I/2 

c(v,  s)  =  (v  vs)'  -  (v  +  Y  s  )  -  vln  Y  - ^ j2 

v+  (v  +  Y  s  )  ' 


(24.1) 
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It  can  easily  be  seen  that 


e(0,  s)  =  s(i  -  -y)  >  0 


e(v,  s)  —  -vln  y  >  0  (as  v  -*•  ®) 


(24.2) 


also 


„  t  .  2  ,  2.1/2 

■^=  -vln v  ~—T~zjrrn  > 0 

v  +  (v  +  Y  s  } 


(24.3) 


Thus,  t(v,s)  is  always  positive  and  increases  monr/conically . 

The  solution  of  equation  (22)  is  now  routine.  As  u,  at  the  back  face 
and  at  0  =  0  is  all  that  is  required,  only  and  will  be  solved  for,  since. 


Ill  _  h(ys)  tt  _  iH  1  T1 
2  ys  1  '  ys  h(qys)  2 


(25) 


Solving  for  and  U^,  one  gets 


Ul  =  - 


f(v,  s)r 


1 


Hs24  h(’1YS) 


4  sinh  6  4v  ,  ,  ,,  ,  .  .  , 

11  h(Y3)'h(YS)  '  “TT?  h(s)h^vs)  sinh  < 


U2  = 


f(v,  s)r 


1  2v 


2.  2  2 

ps  A  y  s 


h(ys) 


4  sinh  6  4v  ,  .  ...  .  .  , 

1  h(YS)h(„s)  '  TT5  h's)h<’lV=)  smh  < 


(26) 
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The  procedure  now  is  to  substitute  equation  (27)  into  equation  (26), 
integrate  on  v  on  an  appropriate  path  to  invert  the  8  transform  and  evaluate 
for  large  values  of  s.  As  it  stands,  this  is  difficult  because  of  the  form  of 
A  and  h,  and  because  the  saddle  path  in  the  v  plane  is  not  yet  defined.  In 
order  to  proceed,  one  recasts  A  into 


c  6  8 


+  “4-5-4  h(s)h(vs)h(r|s)h(riYs)  e"Ce"6j  (28) 

Hs  / 
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or,  retaining  only  powers  up  to  s"  ,  into 


a  - 1  ~  4h(s)h(\s)h(r|s)h(nvs 

t  €  8 

e  e  r| 


j1  +_zirhw(lis)  - 

(  ns  L 


i2(vs)h2(riYs) 


(l  +|l  +  [h2(s)h2(n»)  t  +  e-26) 


+  ~42^  4  ^(s)h(vs)h(r)s)h(r)ys)  e"E  e'6 
*1  Y  s 


Using  equations  (26)  and  (29)  in  (25),  one  gets 


~  f(v,  s)r ,  -  €  -6  /  _  2 

u^-4h(s)h(Ys)— -2-  -  e  ve  fl 

q  p  Y  \  Y  n  s 


ji  +  ^[h2(s)h2(ns)J  + 


1  ,  h2(Ys)h2(Tws) 


X  sinh  6  -  — . h(s)h(Ys)h(ris)h(qYs) 

L  U  Y  s 


(■*!'*£( 


h2(s)h2(r1s)  + 


h_2(Ys)  t  h2^Y»)ll(e-2t  +  e-26j 

Y4  J! 


+  ~4~2~4  h(s)h(Ys)h(Tls)h(riYs)  e"€  e‘5  +  •  •  • 
n  Y  s 


where,  for  repetition, 


2*P0  -s/v  e§(  l»  s/v) 

f(v,s)  =  e  s/vl___ - 

s  h(s/v) 
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In  order  to  invert  equation  (30)  with  respect  to  0,  one  takes  the  inverse 
Fourier  transform 


/CO  . 

5{v,  s)  e  dv  (32) 

This  means  that  integrals  of  the  form 

1=  rP(h)vke-€+g('Vl*s/v)-lv9dv  (33) 

-CO 

must  be  evaluated  for  large  values  of  s.  p(h)  represents  the  polynominals 
in  the  various  h  functions.  The  integrals  are  evaluated  by  finding  the  saddle 
points  of 

-€(v:  s)  +  g(  |v  |,  s/v)  -  iv0  (34) 

However,  in  the  present  case,  interest  is  centered  only  on  the  point  0=0. 
This  simplifies  matters  considerably.  Further,  the  h  function  cannot  be 
simplified  until  the  appropriate  saddle  point  paths  are  determined. 

Some  general  considerations  can  be  inferred  from  equation  (30),  before 
inversion.  The  quantity  c  is  associated  with  dilatational  effects  as  can  be 
seen  by  putting  «  in  a  dimensional  form;  similarly,  6  is  concerned  with 

equivoluminar  effects.  In  order  to  arrive  at  equation  (29),  a  series  of  terms 

c  c  2  €  2  § 

with  higher  powers  of  e  and  e  is  terminated  with  the  term  (e  +  e  ) 

_  c  -§ 

and  e  e  .  These  terms  represent  reverberations.  Further,  in  equation 
(30)  the  terms  e  C  and  e  ^  represent  the  time  delay  necessary  for  the  P  and 
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S  waves  to  reach  the  back  face.  This  becomes  clearer  as  the  analysis  pro¬ 
ceeds.  For  simplicity,  only  the  first  reverberations  are  considered  and 
the  terms  containing  e  6  e  etc.,  are  dropped.  It  is  obvious  how 
to  treat  them.  Since  we  are  only  interested  in  u  directly  behind  the  pressure 
front,  only  the  case  of  0  =  0  will  be  treated.  For  other  values  of  8,  the  saddle 
path  becomes  much  more  complicated. 


IIL  INVERSION 


From  the  above  remarks,  one  can  see  that  during  the  first 
reverberation 


1  f(V,S)rl  h(s)h(Ys) 
u  “  „Z  Ys 

(,  +  4v2 
|  4  4 

r(h2(S)h2hs)  + 

s  n 

v 

X 


A  2 

j-T  h(s)h(Vs)h(ns)h(TiYs)  e€ 

n  Vs 


-e 

e 


-6 

e 


(35) 


where  f(v,  s)  is  given  by  equation  (31). 

Thus,  integrals  of  the  form 

I  =/'"p(h),ke-s/v-'+S<M.  (36) 

where  p(h)  is  a  polynomial  function  of  h,  must  be  evaluated.  First  consider 
£ ,  which,  from  equation  (27.2),  increases  monotonically  with  v .  Next  con¬ 
sider 


*[m-(*)H 

This  function  starts  out  at  v  =  0.  Since  g(v,x)  is  monotone,  decreasing  with 
v,  then  g(v,  s/v)  -  s/v  is  always  negative.  Thus,  the  argument  of  the 
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exponential  in  equation  (36)  is  maximum  at  v  =  0.  A  simple  expansion  shows 
it  has  a  saddle  point  at  v  =  0.  After  some  algebra,  one  has,  for  s  »  v. 


-  g  Mv)M  =  s'1  -  Y>  *£{(?- l)  *  v] 


The  second  terms  in  equations  (37)  and  (38)  are  positive  so  that  u  =  0  is 
always  a  saddle  point.  Since  the  saddle  point  is  at  v  =  0,  the  h(v,  s)  functions 
must  be  expanded  around  v  =  0  for  large  values  ox  s.  By  expanding  and 
collecting  terms,  one  gets 


?  i  (?  ,i/2 

u  -  rrsr  p0  (2'>  ,~i, 


*  *  '■  n2  ‘  nV  +  *)(‘ 1  v‘2 


-6-fg(s/v)-s/s' 

2  2  e 
n  Ys 


-  v|  e-«+g(s/v)-s/v 


In  equation  (39)  terms  of  order  less  than  l/s  are  dropped.  Equation 
(39)  is  now  multiplied  by  dv/2ir  and  integrated  from  -<o  to  +co.  This  yields 
the  velocity  at  8  =  0.  The  result  is 


X  '  2fi  S  [(1/Y  -  1)  -r  v]1/2 


■hi 


4/T12  -  2/V2n2  T  1/4  T  1/4Y2  -  y2/4 
(1/Y-  1)  4  v 


Xe-s(i-V).i  4 

s  n2v 

X  vggL  e-*n(t-V)J  , 

li/n(i/Y  -  i)  +  V]J/Z  ) 

-  -Vjs 

Now,  e  *  ,  in  nondimensional  terms,  represents  the  time  delay 

necessary  for  a  P  wave  to  travel  across  the  specimen.  Similarly,  e 
represents  the  time  for  the  S  wave  to  traverse  the  specimen.  With  these 
times  denoted  by  t^  and  tg,  respectively.  Then,  equation  (40)  inverts  to 


-4rlPo  (yv)1/2 


h 


(x + 2p)  (i/Y  - 1)  j 1  ■  ^ "  yki H(t  ■  y  ■ (t  ■  vk2 H(t  ■  v 


where 


_J_  (i/Y  -lllv 
2  n2v  [i/n (i/Y -  i)  +  v]3/2 


/ 


It  is  thus  seen  that  the  predicted  response,  equation  (41),  has  the 

same  form  as  that  shown  in  Fig.  2;  namely,  a  step  function  followed  by  a 

ramp  associated  with  the  P  wave,  and  then  followed  by  a  downward  ramp 

associated  with  the  S  wave.  An  interesting  feature  of  the  P  wave  ramp  is 

2 

that  it  is  possible  for  it  to  be  negative  or  downward.  If  v  can  be  neglected 
(v  =  0.  1  in  most  practical  cases),  then,  for  an  upward  ramp  to  occur. 


2  2/n  -  1/4  3  -  7v 

Y  2  “9-1 5v 

4/m  t  1/4 


where  v  in  this  instance  is  Poisson's  ratio.  The  physical  meaning  of  this  is 
not  clear.  The  negative  slope  for  low  values  of  Y  could  imply  that  a  spreading 
out  process  of  P  waves  is  occurring.  The  possibility  also  exists  that  the 
approximation  is  poor  in  this  region.  It  is  more  likely  that  a  downward  slope 
is  due  to  the  spreading  out  process.  It  is  worth  recalling  that  Y  =”  1  corre¬ 
sponds  to  a  thin  ring,  while  Y  "=  0  corresponds  to  a  thick  cylinder. 

During  the  next  reverberation,  additional  terms  due  to  the  neglected 
€  c  2^ 

powers  of  e  ,  e  ,  and  e  re  enter  in.  Qualitatively  speaking,  these 
result  in  steeper  slopes  and  additional  delays.  The  actual  analysis  is 
not  gone  into  here,  but  experimental  results  bear  this  out  (Fig.  4). 
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\ 


1.0  /tsec 

Fig.  4.  Multiple  Reverberation  Trace 
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IV.  NUMERICAL  RESULTS  AND  CONCLUSIONS 


Since  this  paper  is  concerned  with  the  results  of  a  particular 
experiment  only  that  case  was  treated  in  the  numerical  studies.  A  para¬ 
metric  study  was  not  carried  out. 

The  specimen  was  made  of  aluminum,  with  an  inner  radius  of  4.  0  in. 
and  a  thickness  of  0.40  in.  It  was  2  in.  in  the  axial  direction  and  2  in.  across 
in  the  circular  direction.  The  material  properties  were  taken  Lorn  hand¬ 
books  as  well  as  inferred  from  the  oscilloscope  traces.  The  results  follow: 

E  =  10.4  x  106  psi 
V  =  0.  33 

p  =  3.  90  X  10^  psi 
w  =  0.  101  lb/in.  3 

c  =  2.45  X  103  in.  /sec 
P 

c  =  1.22  X  103  in.  /sec 
s 

The  maximum  rise  above  the  initial  discontinuous  jump  is  the  quantity 
of  greatest  interest  and  occurs  at  t  =  tg  -  t  .  It  is  calculated  from  equations 
(41)  and  (42)  to  be  a  24.42  percent  rise.  Although  numerical  values  taken 
from  oscilloscope  traces  are  difficult  to  obtain  accurately,  a  sample  of  four 
gives  values  of  20.  7,  20.  5,  21.  6,  and  22.  9  percent.  This  is  an  excellent 


Preceding  page  blank 
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